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We analyze the long time behavior of a quantum computer running a quantum error correction (QEC) code 
in the presence of a correlated environment. Starting from a Hamiltonian formulation of realistic noise models, 
and assuming that QEC is indeed possible, we find formal expressions for the probability of a given syndrome 
history and the associated residual decoherence encoded in the reduced density matrix. Systems with non-zero 
gate times ("long gates") are included in our analysis by using an upper bound on the noise. In order to introduce 
the local error probability for a qubit, we assume that propagation of signals through the environment is slower 
than the QEC period (hypercube assumption). This allows an explicit calculation in the case of a generalized 
spin-boson model and a quantum frustration model. The key result is a dimensional criterion: If the correlations 
decay sufficiently fast, the system evolves toward a stochastic error model for which the threshold theorem of 
fault-tolerant quantum computation has been proven. On the other hand, if the correlations decay slowly, the 
traditional proof of this threshold theorem does not hold. This dimensional criterion bears many similarities to 
criteria that occur in the theory of quantum phase transitions. 
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I. INTRODUCTION 

Quantum computation provides a fundamentally new way 
to process data; as a theory, it is complete and remarkably rich 
fi']. However, any real quantum computer is subject to an im- 
placable physical reality: components of a computer will al- 
ways be faulty due to environmental noise. Hence, the builder 
of a quantum computer faces the conundrum of having to iso- 
late the device from its surroundings and, simultaneously, of 
needing to act on it and read its output |j2|l. Many strategies 
have been devised to address this problem ||li|3LH,|5Llq,|7|l, the 
most general being quantum error correction (ll, [a H Efl fllll • 

Quantum error correction (QEC) should be understood as a 
perturbative approach |,1Z1 . where one can estimate the proba- 
bility of having an "error" in the wave function of the quantum 
computer after a certain time. It is naturally formulated as a 
perturbation expansion in powers of the coupling between the 
computer and the environment lfl2ll . QEC cannot, in general, 
perfectly correct the quantum evolution, and the interference 
of the amplitudes for the various processes that occur implies 
that quantum information is always lost to the environment 
lfT2ll . However, as we discuss below, QEC can very effec- 
tively slow down this loss. In fact, a central theoretical result 
is the "threshold theorem": it states that if the error probability 
is smaller than a critical valu e, q uantum cornputation can be 
sustained indefinitely ||ll[ll[ll[ll[Il[Il[lllMIIltl.The 
word "indefinitely" deserves some clarification: For the prob- 
lems that we discuss, it means that given a calculation and a 
desired precision, it is always possible to construct a quantum 
circuit that will provide the correct result with high enough 
probability. 

QEC has been largely developed using phenomenological 
"error models". Rarely is a connection to a microscopic quan- 
tum dynamica l sy stem found in the literature (see, however, 
Refs. M22i l23i [24l l25ll ). In contrast, here we pursue exactly 
such a connection: We discuss the formal steps needed to link 



the theory of error correction with microscopic Hamiltonian 
models. Furthermore, because of the perturbative nature of 
the method, it is possible to draw a close parallel between the 
"threshold theorem" and the theory of quantum phase tran- 
sitions. We find that if a certain inequality holds, an error 
threshold always exists. When the inequality is not satisfied, 
either a new version of the threshold criterion is required or 
fault tolerant quantum computation is not possible at all. For 
the moment, we are not able to distinguish between these two 
possibilities. 

Our analysis is based on the following assumptions. First 
and foremost, we assume that it is possible to perform the 
building blocks of quantum error correction, namely, prepa- 
ration of states, quantum gates, and measurements. Second, 
we consider that the environment is described by a free field 
theory in which thermal fluctuations can be effectively sup- 
pressed. Finally, the main simplifying assumption of our dis- 
cussion is that the qubits are sufficiently separated in space 
for an entire error correction procedure to be performed be- 
fore correlations between nearby qubits develop. The proba- 
bility of an error in an individual qubit within a QEC cycle is, 
therefore, independent of all other qubits. This does not im- 
ply that there are no spatial correlations; rather, they develop 
on longer time scales, while the error correction procedure is 
done faster than a certain characteristic time. We emphasize 
that this hypothesis is not a limitation of the general theoreti- 
cal framework that we describe, but simply a way to connect 
to the traditional proofs of the "threshold theorem" in terms 
of stochastic error models. 

The paper is organized as follows. Because of the inter- 
disciplinary nature of the subject, this Introduction continues 
with a discussion of two points. First, the difficulties in taking 
into account correlations in the environment are explained in 
Sec. II Al from a perturbative point of view. Then, in Sec. lIBI 
we discuss the QEC method from a physics viewpoint and 
present some results for the standard stochastic error model. 
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We start the body of the paper by developing the relation be- 
tween error models and quantum codes (Sec. The key 
issue of QEC in a correlated environment is treated in Sec. Hill 
Our main results delineating when the perturbative treatment 
is valid appear in Sec. HV] At the end of this Section, we 
provide a brief comparison between our results and those of 
Ref. 26. Sec.|V]discusses parallels between the threshold the- 
orem of QEC and the theory of quantum phase transitions. 
Finally, in Sec.|VT]we summarize our results and comment on 
some open problems. 



A. The problem of correlated environments 

In order to set the stage for the analysis in the presence 
of QEC, we first look at the problem of errors created by 
a correlated environment in an unprotected system. In the 
Schrodinger equation governing the time evolution of a quan- 
tum system, the Hamiltonian H can usually be separated into 
a single-particle term i/o and a many-particle interaction part 
V . A formal solution of this equation is given by the Dyson 
series in the interaction picture. Solution by iteration shows 
that the time evolution operator is 



(1) 



with Tt denoting the time ordering operator and V{t) — 
eh^oty^-j-Hot^ Y represents the interaction between the 
quantum computer and its surroundings, each insertion of V 
in Eq. ([T]i corresponds to an "error" in the computer evolution. 
Hence, Eq. ([T]i provides the natural framework to study the ef- 
fects of the environment on the state of the quantum computer. 

It is always possible to give an upper bound to the "error 
probability" [27]. The reason is that Dyson's series is abso- 
lutely convergent for finite times and bounded operators (see 
Appendix IaIi. In short, the bounding is done by defining the 
"sup" operator norm and the evolution operator with at least 
one "error" (one insertion of V), 



£{t) = U{t,0) - 1 



dt'V{t') Uit',0) . 



(2) 



The norm of £ is related to the probability of having errors in 
the computer. The calculation is simple and yields 



\£im< 



1 



dn\v{n\\< 



At 



(3) 



where we used the triangular inequality, the unitarity of U, and 
defined A as the largest eigenvalue of V (with corresponding 
eigenvector ^'a). One can understand this bound as simply a 
restatement of Isinxl < \x\, as follows: 



1 /■* 

I-Ttcos- / dt'V{t') 



(4) 




1^ < • (5) 



The norm \ \£\ \ has been very useful in problems involving 
non-Markovian noise ll26l IztI l28l l29i i3Q] . However, in QEC, 
an analysis based on the bound Eq. Q only makes sense when 
||£|| ^ 1, while we are concerned with the long time limit, 
I At I ^ 1, for which this bound on the norm of the error di- 
verges. In this case, Dyson's series is only asymptotically con- 
vergent and the "sup" norm is of no practical use. Hence, it is 
important to express the error probability differently. 

We must go back full circle and reexamine the Dyson series 
for the time evolution of a particular state, instead of the worst 
case scenario explored by the "sup" norm approach. Hence- 
forth, we will be mainly interested in an interaction Hamilto- 
nian with the general form 



(6) 



where A ^ 1 is a coupling constant, L is the size of the sys- 
tem, and / is some function of the degrees of freedom of a 
free theory whose Hamiltonian is Hq. Because we are inter- 
ested in correlated non-Markovian noise, we assume that the 
free fields are such that the asymptotic expression for the two- 
point correlation function is a power law. 



(Vl/|/(xi,ti)/(x2,i2)|*)^^ 



1 



1 



_ (7) 

where Ax — |xi — X2I and At ~ \ti — t9| ll3lll . Here, 6 is the 
scaling dimension of /, z is the so-called dynamical exponent, 
and I is a fixed eigenstate of Hq (which we will usually take 
to be the ground state of the environment). 

The motivation for developing a perturbative expansion of 
the evolution operator (the Dyson series in the interaction pic- 
ture) is the hope that a few terms in the series or a summable 
family of them will capture most of the physics. It is then 
assumed that small coupling can guarantee fast convergence. 
However, since is not necessarily small, the number of 
terms that contribute substantially to the series can grow faster 
that the smallness of consecutive terms. In order to see that, let 
us calculate the probability of an evolution with errors using 
Eq. dUl, 



(vl/|£t it)£it)\^) 



= 1 



(*|T^cos 



1 



- / dt'V{t') 



(8) 

Since we are assuming a non-interacting free Hamiltonian, we 
can use Wick's theorem. It is then straightforward to show 
that there is at least one term at each order m in the series that 
contributes "extensively" as ^ A^'"(Lt)^'"(^+^~'^). A simple 
example is given by the series of "bubble" diagrams, where 
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the m order term is given by the contractions 

/ dh... ' dt^ {Vih)V{t2)) ... {V{t,n-l)Vit„,)) . 

Jo Jo 

(9) 

Disregarding numerical prefactors unimportant for our discus- 
sion, we sum the series as a geometric progression to obtain 

2 ^ 1 + A2 {LtyiD+.-S) ■ ^ > 

Therefore, for D + z — 5 > Q there is no guarantee that the 
perturbation series converges. Conversely, if Z) + z — 5 < 0, 
higher-order terms in the series should be increasingly less 
important. Thus, for D + z — 5 > Q the probability of an evo- 
lution with "errors" tends to one, whereas for D + z ~ 5 < 
it will depend only on the "non-extensive" terms in the se- 
ries. The same analysis can be immediately transported to the 
study of the fidelity \{'^\ U{t) |^)|, where we see that for a 
relevant perturbation, D + z — 5 > 0, the overlap between 
the initial state and the evolving wave function tends to zero 
(an orthogonality catastrophe). This sort of "infrared" prob- 
lem provides a contact point with the theory of quantum phase 
transitions, where the same kind of considerations also appear 
when calculating the partition function using the imaginary 
time formalism (see AppendixlBl). 

In the body of this paper, our main goal is to transfer these 
ideas of relevance and irrelevance of a perturbation to the evo- 
lution of a quantum computer protected by QEC. 

B. Quantum error correction 

Quantum error correction is arguably the most versatile 
method to protect quantum information from decoherence 
ll32ll . It is a clever use of two features of quantum mechanics: 
entanglement and (in its traditional form) wave packet reduc- 
tion due to measurement. Thus, before we start our discussion 
of QEC, it is important to carefully define what we mean by 
entanglement and decoherence. 

An entangled state of two quantum systems is a state that 
cannot be described as a direct tensor product of states of in- 
dividual systems or probabilistic mixtures of tensor-product 
states. As an example, consider two physical qubits (hereafter 
referred to by the subscripts 1 and 2). Each qubit has a Hilbert 
space isomorphic to a complex projective plane of dimension 
one, CP^ (see Appendix ICl for details). However, the com- 
bined Hilbert space is not isomorphic to CPj^-) x CP(2), but to 

the much larger CP^. All states in CP^ outside CPj^j x CPjj) 
are said to be entangled. An important subtlety is the implicit 
notion of a preferred "basis". Although we can choose from 
an infinite number of CP^ x CP^ subspaces inside the same 
CP"^, nature gives us a natural choice, namely, CPjj^-) x CP(2) ■ 
In the working of a quantum computer, entanglement has 
two opposite roles. On the one hand, entanglement between 
qubits is the key element in a quantum cornputation that dis- 
tinguishes it from its classical counterpart II33|] . On the other 
hand, when the computer and the environment become en- 
tangled, precious quantum information is lost. Usually, the 



latter effect is referred to as decoherence. In the literature, 
there are two different definitions of decoherence. In a strict 
sense, decoherence is the decay in time of the coherences (off- 
diagonal elements of the reduced density matrix), while dissi- 
pation involves the exchange of energy with the environment 
and changes in populations (the diagonal terms of the den- 
sity matrix). However, the word "decoherence" is also used 
in a broader sense involving changes in both diagonal and 
off-diagonal entries of the density matrix. In this paper we 
choose the latter use of the word. The reason is that from a 
quantum error correction perspective changes in diagonal and 
off-diagonal entries are "dual" to each other IJJ]. 

There is a simple heuristic explanation for error correction: 
Usually, noise is regarded as a local phenomenon, thus its 
damaging effect in the computer should be less pronounced 
if the information is delocalized among several qubits. This is 
precisely how classical error correction codes work. A simple 
example of the latter is a majority vote, where the information 
of a bit is copied into three physical bits, 0^000 and 1^111. 
If the probability of an error in a given qubit is e, the probabil- 
ity of having two independent errors, and consequently a total 
information loss, is ^ e. Thus encoding increases the level 
of protection of the information. 

It is tempting to start explaining QEC from this perspec- 
tive. However, the no-cloning theorem ^ states that it is im- 
possible to copy an unknown quantum state. The alternative 
approach is to use an entangled state involving two or more 
qubits to store the quantum information. This clearly delocal- 
izes the information, but it is at odds with the intuitive notion 
that entangled states are in general more fragile to the effects 
of the environment (this intuition is driven by the quantum-to- 
classical transition due to decoherence, see Appendix iDl for a 
concrete example). Thus, delocalizing the information using 
entanglement does not alone solve the problem. It is possible 
to use unitary operations to transfer the entanglement between 
the qubits and the environment to a constant fresh supply of 
ancilla qubits (T,^3^. However, it is more traditional in QEC 
to use the partial measurements of some ancilla qubits to re- 
duce the quantum interference with the environment iH] . Mea- 
surements here have to be understood as the projection of the 
state of one of the qubits (an ancilla) onto a certain basis or 
reference state. The outcome of this projection is a classical 
bit ("zero" or "one") and is called a syndrome. The partial 
wave packet reductions caused by syndrome extraction steer 
the long-time evolution of the quantum computer Recently, 
it has been shown that the duration of the measurement is not 
fundamental to the QEC procedure jSSt]. In fact, this process 
can be quite long without jeopardizing the method. 

A simple example illustrates how QEC works lU |^ [T^ . 
Suppose that we have an error model consisting of indepen- 
dent baths for each qubit which can cause only phase errors, 
and an initial qubit in the state \iIjq) — a\1) + f3\i) that we 
want to protect. The 3-qubit code provides the simplest error 
correction procedure for this problem. In Fig. [T] we define the 
encoding/decoding methods in a QEC cycle. At the end of a 
cycle, the probability of measuring the syndrome of a phase 
flip error in one of the three physical qubits is L39il 

Pi = 3e, (11) 
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FIG. 1: A3 qubit quantum error correction (QEC) code Ql, |^ 
HE m, El- The initial wave function, l^o) (|T> + |i>)/2 (g) 
(IT) + li>)/2, is encoded by two controUed-NOT (CNOT) gates, 
-RcNOT = a~'afaj + af , into an entangled state IV'encode) = 
a It) + /3 ID with ]f ) = (1TTT> + ITii> + lin> + liiT)) /2 and 
ID = (liU> + UTT) + lnT) + ITTJ.))/2. After some time, it is 
decoded by a second pair of CNOT gates. An error in \ip) is iden- 
tified by measuring the values of erf and erf (rectangle). The QEC 
cycle ends with the correction of a possible phase-flip (arrow). 

and the probability of the syndrome indicating no error in the 
logical qubit is 

Po = l-Pi- (12) 

The residual decoherence that can not be corrected by the 
QEC procedure is closely related to these probabilities. In the 
case of a cycle in which the syndrome indicates that one error 
occurred in any of the physical qubits, dephasing of the logi- 
cal qubit is give n by the reduction of the off-diagonal density 
matrix element ll39tl . 

« a/3* (1 - 2e) , (13) 

while for a cycle with a syndrome indicating no error, the de- 
phasing is weaker, 

pfl^al3*{l-2e^). (14) 

After N of these cycles, the probability of having m uncorre- 
lated errors is 

with an associated residual decoherence of 

(l-2e3)^-"(l-2e)'". (16) 

An elegant visualization of these events is given by a "syn- 
drome history diagram" of Fig. |2] (see for instance Ref . |28| for 
a similar discussion). An ordered set of syndromes labels a 
particular evolution of the logical qubit. From the syndrome 
history one can find the most likely evolution and the asso- 
ciated residual decoherence. For our 3-qubit code example, 
the most likely evolution is given by the mean value of m, 
in = Npi. Thus, the residual decoherence of the logical qubit 
is given by 

p^i «a/3*e-6^^'. (17) 

Therefore, as long the number of QEC cycles N ^ e^'^, the 
probability of measuring the correct initial state of the logi- 
cal qubit is very high. We can quantify the amount of infor- 
mation that is lost by calculating the von Neumann entropy 
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FIG. 2: A syndrome history diagram. Each solid line represents the 
evolution of a logical qubit. At the end of a QEC cycle, a phase flip 
error is detected or not with probabilities pi and po, respectively. A 
path provides the history of the logical qubit and is recorded as a 
sequence of syndromes. 



S = -tr(plnp): 

lim S « l2N\a\^\P\\^\l-lii(l2N\a\^\l3\''e^)](l8) 
lim S « -\afh-i\af-\f3f\n\f3f. (19) 

Note that the loss of information can be substantial if the num- 
ber of cycles is so large that N ^ e^^. 

If the information needs to be protected for a long period 
of time, we have to modify the protection scheme. The most 
straightforward approach is to consider a concatenated circuit 
where each qubit in Fig. [T]is a logical qubit itself and each 
gate is a logical gate, resulting in an effective reduction of pi. 
Layers and layers of protection can be added as needed 1 1, 29]. 
A chief concern when applying this approach is whether the 
steps required in the addition of more qubits and operations do 
not actually increase the chance of errors (since they increase 
the combinatorial factors in the probability distribution). This 
question is addressed by fault-tolerant quantum computation 
theory ill [H [M [HIS 29], which has as its main result 
the so-called threshold theorem: If the "noise strength" e is 
smaller than a certain critical value, then the introduction of 
an additional layer of concatenation improves the protection 
of the information. 

A key ingredient in the derivation of the noise threshold 
is the assumption that a probabilistic structure similar to the 
one that we outlined above exists. Here rests the main con- 
cern of this paper There are many physical situations where 
an environment can induce strong memory effects and spatial 
correlations among qubits. Hence, it may not be obvious how 
to define the "error probabilities" of a qubit. This hinders the 
traditional theory of QEC and threshold analysis, thus moti- 
vating a careful study of the dynamics of quantum computers 
protected by QEC. 
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II. ERROR MODELS AND QUANTUM CODES 

The syndrome history used to describe the logical qubit his- 
tory can be converted into a more formal description of the 
computer dynamics. In our discussion, we will assume an en- 
vironment. Ho, described by a free field theory with an ultra- 
violet cutoff A, a characteristic wave velocity v, and a dynami- 
cal exponent z. Although simple, a free field theory faithfully 
represents many physically relevant environments: the elec- 
tromagnetic field, a phonon field, spin waves, a bosonic bath, 
or, more generally, any two-body direct interactions between 
qubits that was split by a Hubbard-Stratanovich field. In addi- 
tion, we include in the Hamiltonian a term to account for the 
sequence of quantum gates performed on the qubits, HQc{t). 
Hence, the total Hamiltonian is 



H{t) - Ho + Hcic{t) + V. 



(20) 



The interaction term will be assumed to have the form of a 
vector coupling between qubits and the environment. 



A. 



(21) 



a={x,y,z} 



where (f{x) are Pauli matrices for the qubit located at x, 
are the coupling strengths, and fa (x) are functions of the en- 
vironment operators [38]. Since [Hq, Hqc] = 0, we adopt an 
interaction picture that follows not only the environment but 
also the evolution of the computer (see Appendix |E|i. In this 
rotating frame, the evolution operator is 



C/(i,0) ^Tte 



-ilodt'vit') 



(22) 



The interaction V{t) depends on the quantum code and its 
implementation. Nevertheless, there are two possible ways to 
keep the discussion code independent: 

(i) In our previous work [39, 40], we assumed that quantum 
gates were performed faster than the environment response 
time (which is of order the inverse of the ultraviolet cutoff 
frequency A). We call this approximation the "fast gate" limit. 
For this case, we have the evolution of the computer between 
gates given by 



V{t) = 



E 



E 



^/„(x,t)a"(x), 



(23) 



} 



with /q,(x, = eTi^°*' fa{'x)e~^^°* . Then, when a gate is 
performed the action on the qubit is instantaneous and the sub- 
sequent evolution is once again governed by Eq. ( |23] ). 

(ii) A second possibility is to derive an upper bound on the 
effects of correlations. In order to do that, we must first dis- 
cuss how slow gates, which are performed over time intervals 
larger than Tc ~ 1 /A, change Eq. (l23T l. Then, we can define an 
effective interaction Vcs that takes into account the slowness 
of the gates and serves as an upper bound to the exact (and 
code-dependent) V . Clearly, the real experimental situation 
rests between the two limits (i) and (ii). 

Before we begin a detailed description of how to handle 
case (ii), let us note that here the terminology "fast" and 



"slow" gates follows the QEC literature: Fast (slow) gates 
have a duration much shorter (longer) than t^. However, as 
will become clear later, the relevant time scale that appears in 
the study of correlation effects is the period or duration of the 
error correction cycle, A. Thus, in that context, short ("fast") 
or long ("slow") dynamical effects will be naturally defined 
with respect to A, and not to Tc- 

Any quantum computer code is just a rotation in the Hilbert 
space of the qubits and can be described as a trajectory 
on CP^^~^, where N is the total number of qubits. In 
the Schrodinger picture, the evolution is given by the natu- 
ral action on S"^^^^ by SU{2^). The most general fault- 
tolerant quantum circuit is therefore defined by the Hamilto- 
nian HQc{t) = J2 ^i(^) Cj", where {ej} are the generators of 
the Lie algebra of SU{2'^). The evolution operator associated 
with this Hamiltonian satisfies the integral equation 



Wit,0) 



1 



ft JO 



dt'HQc{t')W{t',0) 



(24) 



such that the computer state vector at time t is given by 
\ip(t)) = W{t,0) IV'(O)}, where |V^(0)) represents the initial 
state of the computer Therefore, in the interaction picture, the 
interaction operator is given by 

V{t) = W''it)ei""^Ve-T:""^W{t) (25) 

= E E yI"*"'"^-^''^"'". 

X a—{x,y,z} 
^ ot.—{x.,y,z} 

Since W{t) is a SU{2^) matrix, then 

G" (x, t) = (t) a" {^)W{t) (27) 
is another matrix of SU (2^), and we can write 

^w-E E Y^"^^'*)^"^'''*)- ^^^^ 

X a — {x,y,z} 

Although the expression in Eq. ( |28] l is general, it is not very 
instructive. Furthermore, it is very undesirable from an er- 
ror correction standpoint: since G{t) is an arbitrary matrix of 
SU{2^), the V{t) in Eq. ( |28T l in principle generates a highly 
complex correlated error that is nevertheless first order in the 
coupling to the environment. The problem with the deriva- 
tion of Eq. (l28T l is that it is too general since we assumed that 
arbitrary rotations are performed at each single step. How- 
ever, one of the cornerstones of quantum computation is that 
such general rotations can be approximately decomposed into 
a series of elementary gates [1]. Hence, our strategy will be 
to specialize the calculation to these elementary gates and as- 
sume that general rotations can be implemented by a finite 
series of such gates which are well resolved in time. 
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A. Single-qubit operations 

When only single-qubit operations are performed, we have 
HQc{t)^Yl E ^a(x,i)a"(x). (29) 

In this case, W{t) is the product of SU{2) matrices acting in 
each qubit's Hilbert space. Thus, G"(x, t) simplifies to 



G?(x,t) 



Pie 



— i(p 



a"(x) 



Pie 



(30) 

where Pi + pi ^ 1 ™d {Pi' P2i '/'i ^6 functions of x and 
t. The single-qubit rotations yield an expression of the form 

G?(x,i)= J2 5"''(x,t)a^(x) (31) 

/3={l,x,y,z} 

for some g"^{x,t). By decomposing the operators fa and 
functions g"'^ into their Fourier components, we can give a 
more formal meaning to "fast" and "slow" gates, 

/„(x,i)g"'5(x,i)= e^^"^+"^^Va(x,u;i)g"''(x,c.2). 

a)l|<A.[^2 

(32) 

Hence, if we define = cji + ti>2i we can rewrite the pertur- 
bation as 




a/3 , 



Cr/3. 



(33) 

In the limit of fast gates, \uj2\ > A, / and g are not convolved, 
since they have distinct frequency domains. Therefore, the 
noise operators fa are unaltered by the rotation. However, if 
g has a significant weight at frequencies smaller than A (slow 
gates), one must convolve / with g, yielding a substantially 
different noise operator. 



B. Two-qubit operations 

The general Hamiltonian for two-qubit gates is of the form 

x,y a,p={x,y,z} 

(34) 

However, one can also generate a full set of gates using instead 
a single type of interaction, 

where a and b are fixed for each gate (x, y). In order to see 
that this is sufficient we can for instance set a = 6 = 0. This 
generates the liquid NMR Hamiltonian |4ll], where the Ising 
interaction, Eq. (|35] |. and single qubit rotations can be used to 
generate a control gate. 



We keep a and b arbitrary. However, for the sake of sim- 
plicity, we assume that only operations between disjoint pairs 
are allowed; that is, if J (x, yi,t) ^ 0, then J (x, y2,t) = 
for all y2 ^ yi- It is then straightforward to write down W{t) 
in a compact form: The time-ordering [Eq. (I24l i1 is automati- 
cally taken care of by the sequence of gates, while for a gate 
involving qubits x and y the contribution to W{t) is 

W{x,y,t) = cos[0(x,y,t)] (36) 
+ ism [0(x,y,i)]a'^(x)a''(y), 

where 0(x, y, t) = dt'J{x, y, t'). Hence, a two-qubit rota- 
tion yields 

G^(x,t) - sin[20(x,y,i)]6'^"V^(x)a''(y) 

+ cos[2(?(x,y,t)](l-,5a,„)a"(x) 

+ '5a,„a"(x), (37) 

where e'^"^ is the usual antisymmetric tensor. 

The first term on the r.h.s. of Eq. (IJTT i tells us that the 2-qubit 
gate can propagate the error from the qubit at x to the qubit at 
position y. However, it also tells us that it is possible to choose 
a particular gate where this propagation does not happen (by 
choosing a = a, for instance). Unfortunately, propagating 
errors in the quantum circuit is in general unavoidable (since 
the only gate that commutes with all Pauli operators is the 
identity). 

The second and third terms on the r.h.s. of Eq. ( l37T i are 
much less dramatic. They simply describe a local noise that is 
not propagated by the gate. 



C. Upper-bounds for the evolution 

In Eqs. dJTT l and ( |37] |. we showed that one- and two-qubit 
gates can introduce what is seemingly a very complicated 
noise structure. The expressions depend on how the gates are 
implemented, thus hiding a general assessment. We can ad- 
vance the discussion by recalling that W is always an unitary 
matrix. Hence, the coefficients in Eqs.dJTIi and ( |37| i have mod- 
ulus equal or smaller than unity. A suitable upper bound on 
the effects of slow gates is then provided by setting all these 
coefficients equal to one. Thus, the operators expressed in 
Eqs. ( I3TI 1 and ( |37] | gain the upper bounds 

G?(x) - J2 '^"W' ^38) 

fj={x,y,z} 

G^(x) = a"(x) +e'^"''a'^(x)cr^(y). (39) 

G2 still looks troublesome, since it tells us that an error 
in qubit x is propagated to y. However, this is not a prob- 
lem of the finite gate time operation, since an instantaneous 
and perfect gate will also propagate the error in a similar fash- 
ion. In order to obtain an upper bound for the effects intro- 
duced by the two-qubit gates, we precisely follow this fact. 
We consider that all the qubit components are exposed to all 
the noise channels all the time. Thus, we replace Eq. ( |39] l by 
G2 (x) ~ J2i3={x y z} "''^(x) and assume that two-qubit gates 
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are performed instantaneously. In summary, we reduce the 
problem of finite time operation of the two-qubit gate to the 
problem of a noisier qubit environment and propagating errors 
in the quantum code by perfect gates. Now we can rely on the 
theory of fault-tolerance [1, 29J, and simply assume that the 
error propagation is handled by the quantum code. 

The final conclusion is that an upper bound estimate on the 
effects of slow gates is obtained by the interaction Hamilto- 
nian 



E 



E 

={2;, a,. 



A 



/eff(x,i)a"(x,i), (40) 



where 



/cff (x, = T2 



I3={x,y,z} 



(41) 



and A — y Xl/3={a: y z} '■^^ coupling parameter 
Although this is a brutal approximation, it will be sufficient 
for our discussion. As we will argue later, for the purpose 
of determining the effect of long-wavelength correlations on 
the threshold theorem, the only relevant aspect of the is 
their scaling dimension. Since dim /os is in general equal to 
min (dim it is sufficient to use Eq. (|40] l as the worst case 
scenario. 

Thus, in both limiting cases, fast and slow gates, we ar- 
rive at the same functional form for the effective interaction. 
Hence, both cases can be handled simultaneously, and we pro- 
ceed to the analysis of QEC in the presence of this interaction. 
In order to simplify the notation, we hereafter drop the sub- 
script "eff" from the slow-gate operators. 



One of our goals is to find out when it is appropriate to safely 
neglect such uncorrectable errors in the presence of correlated 
environments. 

In hindsight, it is not hard to understand the benefits of 
QEC. Thus, for the sake of readability, we present first a qual- 
itative argument that captures the overall discussion. 

As we defined in the introduction, there are two quanti- 
ties that we are interested in calculating: (i) the probability 
of a given evolution, and (ii) the reduced density matrix of the 
computer. Both quantities are written as a double series in the 
coupling with the environment. On the one hand, the initial 
ket of computer and the environment, evolves in the time 
interval [0,<] by the time ordered series U{t). On the other 
hand, the bra evolves in time with the anti- time-ordered 
series It is only a subset of each series that enters in the 
evaluation of either the probability or the reduced density ma- 
trix, because of the measurements present in the traditional 
formulation of QEC. Hence, it is usually a non-trivial task to 
calculate the necessary expectation values. 

Because we are dealing with a double series, it is natural 
to use a formalism analogous to a time-loop expansion ll42ll . 
There are six (interrelated) Green functions in such an ex- 
pansion: The usual advanced and retarded functions for the 
time-ordered series; the advanced and retarded functions for 
the anti-time-ordered series; and the lesser and greater func- 
tions, which contract a term from the time-ordered series with 
another one from the anti-time-ordered series. This formal- 
ism is often referred to as the Schwinger-Keldysh approach 
||43[ I44I1 . It is usually represented graphically by a double 
contour in time (see Fig.|3]l. The upper leg stands for the time- 
ordered evolution for the time interval [0, i], while the lower 
leg stands for the anti-time-ordered evolution in the reversed 
interval [t, 0]. 

Let us for the moment assume that a short-time expansion is 



III. QUANTUM ERROR CORRECTION IN 
CORRELATED ENVIRONMENTS 

A QEC code is defined as the combination of encoding, 
decoding, and recovery operations. Since we were able to 
make our analysis code independent, the unitary component 
of the QEC protocol is described by C/(A, 0), Eq. (|22]), with 
the appropriate V{t) discussed in Sec. The final ingredient 
in standard QEC is just the syndrome extraction P, which is a 
projective measurement. 

In Ref. 39 it was demonstrated how to define P and its ef- 
fects on U for stabilizer error correction codes. It is important 
to remark that an error which keeps the computer in the logi- 
cal Hilbert space can never be corrected by QEC. This is sim- 
ply a statement that for the general assumptions we make, the 
problem of protecting quantum information never satisfies the 
second criterion of Lafflame-Knill [12] for perfect QEC. In 
simple terms, the criteria states that all allowed errors must al- 
ways take the logical one and logical zero to orthogonal states 
[Eq. (20) of Ref. 12]. By construction, these errors are high- 
order events in the coupling with the environment. Neverthe- 
less, as we already know (II Al l, this fact per se is not enough 
to ensure that such errors will not be relevant at long times. 




(c) (d) 



FIG. 3: Graphical representation of a few fourth-order terms in a 
"time-loop" expansion for either the probability of a given evolution 
or the reduced density matrix (spatial dimensions are suppressed for 
clarity). Points of interaction with the bath (circles) are connected by 
propagation of the environmental modes (wiggly lines). In the top 
diagrams, the time integrals are unconstrained, as would be the case 
for unitary evolution. In the bottom diagrams, the detection of an 
error by a QEC protocol forces the interactions with the bath to occur 
at the same times on both the forward and backward legs in order 
that U and U'' correspond to the same syndrome. This additional 
constraint introduced by QEC is crucial in the long-time behavior. 
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valid and focus on a single qubit. Then, the evolution operator 



In Fig.[3]we represent graphically a few terms of order A"'. All 
of these terms are the product of a second-order term from Ui 
and a second-order term from u\ [see Eq. (l42Ti1. Hence, they 
correspond to two "errors" in the qubit evolution and involve 
the expectation value 

(*! /t (x, i)4(x, <')/a(x, i")/Mx, n Ivf) . (43) 

Using Wick's theorem, we can immediately write (|43] | as a 
product of the non-interacting Green functions. Each possible 
set of contractions leads to the different "diagrams" in Fig. [3] 
We usually do not know when an "error" occurs; hence, 
each Green function is accompanied in the series by a double 
integral in time. This is precisely the case in an unprotected 
computer's evolution or inside a QEC cycle [see Figs|3](a) and 
(b)]. However, a dramatic change happens in a Green function 
between terms for different cycles. When the syndrome shows 
that a particular error occurred in a certain QEC cycle, we can 
re-write Eq. (l43T l to reflect this knowledge: 

(*! /t (x, i)/t(x, t')/a(x, t + 5t)fp{^, t' + 5t') I*) , (44) 

where 5t and 5t' are time variables with range smaller than 
the QEC period. After integrating the "high frequency" part 
(the 5t and 5t' variables), we end up reducing Eq. (|44] | to 

(vI/|/t(x,t)4(x,t')/a(x,t)/Mx,i')l*) (45) 

with t and t' representing a coarse-grained time scale of order 
the QEC period [see Figs. [3](c) and (d)]. Therefore, although 
we are considering terms of the same order in A, the number 
of "time integrals" in the coarse-grained scale (low frequen- 
cies) is half that in the original microscopic calculation (high 
frequencies). 

The simple dimensional analysis of Sec. II Al tells us now 
that QEC has changed the criteria for the stability of the per- 
turbation series at long times. As we demonstrate now, it is 
less stringent than the naive expectation. 

A. Quantum evolution steered by QEC 

It is reasonable to assume that at the beginning of the com- 
putation the computer's state vector, ^/'o, and the environ- 
ment's, (po, are not entangled, 

I* [t = 0)) - l^o) ® \vo) ■ (46) 



for that particular qubit within a QEC cycle is given by 



(42) 

I 

In a realistic situation, ipo would have some initialization error 
and be entangled with the environment to some degree (both 
of which would yield errors in ijj). However, here we neglect 
these effects in order to keep the discussion focused. 

Just as in the case of the 3-qubit code, by the end of a QEC 
cycle the computer will have evolved according to the unitary 
operator C/(A, 0). Then, the syndrome is extracted and the 
computer wave function is projected, 

P„{/(A,0)|*(0)), (47) 

where m corresponds to a particular syndrome, with 
= I and — Pjn- In the case of many logical qubits 
evolving together, then m denotes the set of all the syndromes 
extracted at time A. The last step in the code is the appro- 
priate recovery operation, Rm, depending on the syndrome 
outcome, 

l^'(A)) =i?„,(A + 5,A) P„,[/(A,0)|*(0)). (48) 

Since in a fault-tolerant error correction scheme the infor- 
mation is never decoded (in contrast to the 3-qubit code dis- 
cussed above), the quantum information always remains pro- 
tected. Therefore, we can deal with our two limiting cases 
(slow- and fast-gates) in two different ways. In the case of a 
slow-gate recovery, we formally include it as the initial step 
of the next QEC period. Conversely, in the case of fast gates, 
we assume that the recovery is performed flawlessly in a very 
short time scale after the projection. For the sake of clarity, 
we choose the latter below. We emphasize that this does not 
restrict our discussion, since it is known that the time of re- 
covery is irrelevant to the error correction. In fact, it can be 
postponed all the way to the end of the calculation [35.1 . 



B. Probability of a syndrome history and the loss of 
information 

The first quantity to discuss is the probability of measuring 
a particular syndrome at the end of the first QEC step, 

Vrn - (*(0)| C/t(A, 0) Prn (/(A, 0) |*(0)) . (49) 

The corresponding reduced density matrix is 



a={x,y,z} 

\,Cdtf dt' J2 

a = {x,y,z} 



Jo 



Uk, t)fp{^, t') a"(x, t)a^{^, t') + O(A^) 
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tre [(r-1P„t/(A,0)|^(0))(vl/(0)|j7t(A,0)P„|s)] _ (^o| [{jjo] A, 0) P,„, \s) {f\ P,^U{ A, 0)\^o)] |yo) 



(*(0)|{/t(A,0)P,„C/(A,0)|^'(0)) 



where r and s denote states in the computer Hilbert space and 
tr^ is the trace over the environment Hilbert space. It is possi- 
ble to quantify how much information was leaked to the envi- 
ronment by calculating the von Neumann entropy 

5(A) = -tre [p™(A)ln|p™(A)|], (51) 

where trc is the trace over the computer Hilbert space. 

In Eqs. ( |49l ) and ( |50l ). one clearly sees the important role 
played by the projection operators in the quantum evolution 
steered by QEC. The careful construction of the encoded 
states combined with the measurement (syndromes) reduces 
the quantum interference between different history paths of 
the computer. By partially collapsing the wave function of the 
computer, this traditional form of QEC reduces decoherence. 

Equations i49[ and dSOl ) define the local components of the 
noise. When spatial correlation between qubits can be ig- 
nored, they are related to the stochastic probabilities and den- 
sity matrix discussed in Sec. II Bl [see Eqs. ( fTSl ) and (fT6]l1. 

The generalization to a sequence of QEC cycles is straight- 
forward iH, 

Tw - w»„(iVA,(iV- l)A)...z;^,(A,0), (52) 

where w is the particular history of syndromes for all the 
qubits and 

v^^{jA,ij^l)A) = 

R^^{j{A + 6),jA)P^^U{jA,{j - 1)A), (53) 

is the QEC evolution after each cycle. Each history comes 
with the associated probability 

V (Tw) = (^ol (V^ol I^o) I^o) . (54) 

Finally, there is always some residual decoherence which can 
be found from the reduced density matrix 

, (^o|[(^o|Tt,k-)(r1Tw|V>o)] \^o) 

Pr,s (Tw) = 1 , (55) 

{ipo \ {ipo\ T^Tw 1-00) |</5o) 

with r and s being elements of the logical subspace. This in 
turn yields the entropy 

S (Tw) = -tr, [p (Tw) In \p (Tw) | ] . (56) 

In the following, we will show for Eqs. ( |54] | and (l55b how 
to separate the effect of correlations between different QEC 
cycles from the contributions due to the local component of 
the noise, as defined by Eqs. ( |49] l and ( fSOl l. 



((^ol (^ol ;7t(A,0)P™[/(A,0) l^o) |<^o) 

, (50) 

I 

IV. PERTURBATION THEORY AND THE 
HYPERCUBE ASSUMPTION 

There is one additional issue that we must deal with before 
we can move forward. In principle, even the first order term 
in Eq. ( l42b is akeady beyond the QEC approach that has been 
outlined so far. The reason is that when calculating P or p 
we generate pair contractions of the type (/q(x, t)fa{y., t')). 
Therefore, the probability of finding an error at a given qubit 
is conditional on what happens with all other qubits. This 
automatically hinders the simple probabilistic interpretation 
of QEC that we used in Sec. US 

The fact that we do not want to deal with such conditional 
probabilities leads us to the single most important simplifying 
hypothesis of our work: We assume that the qubits are sepa- 
rated by a minimum distance 

e=(vA)l/^ (57) 

where v is the excitation velocity and z is the dynamical ex- 
ponent of the theory describing the environment. Hence, for 
all X ^ y and \t-t'\<A, we have (/„(x, t)/„(y, t')) « 0. 
It is then possible to assign a probability for the short-time 
evolution of each qubit independently of all others. 

To further organize the analysis we order the qubits in a D- 
dimensional array that defines hypercubes of volume A x 
(see Fig. nil. In summary, for times smaller than A, each qubit 
has a dynamics independent from the other qubits, hence re- 
sembling a quantum impurity problem. However, for time 
scales larger than A, spatial correlations among them are 
present, thus making the problem similar to a spin lattice. 

Ideally, we would like to decompose the evolution operator 
in inter- and intra-hypercube components, 

C/(A,0) = C/<(A,0)C/>(A,0), (58) 

where < labels frequencies smaller than A^^ and > frequen- 
cies in the interval [A^^, A] . Whenever this is possible, we 
can integrate the intra-hypercube part in order to define a "lo- 
cal" evolution and, consequently, a local error probability. 




FIG. 4: Two neighboring hypercubes in space-time, each one con- 
taining a qubit. 
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There are simple noise models where this can be done exactly 
ll39ll . however, in general, this separation is only possible in 
a perturbative expansion. Keeping just a few terms in pertur- 
bation theory is not always adequate, and we must try to find 
ways to improve it. 



A. Perturbation theory improved by RG 

Our objective in this section is to define an effective evo- 
lution operator that can reasonably describe the evolution of 
the qubit within each hypercube. All terms consistent with the 
same syndrome and having the same leading long-time prop- 



erties should be included. Within a hypercube, the environ- 
ment induces interaction of a qubit only with itself; commu- 
nication between qubits at longer times is treated in the next 
section. 

We use the renormalization group (RG) ll45ll to sum the 
most relevant families of terms in the perturbation series. In 
order to improve the lowest order terms in the perturbation 
theory through RG, we need to introduce the next higher-order 
terms in the perturbation series. However, as we discussed 
previously, we are not interested in the full unitary evolution, 
but rather the projected terms obtained after the extraction of 
the syndrome. Therefore, in order to apply RG to the first- 
order term, we need to consider 



J 



UQ(xi,Aa) « ~2h'^°' / dtfa{-Xi,t)- — \eai3'y\Xt3X'r(Ta{A)Tt / dti dt2 //3 (xi,ti) (xi,t2) 0-/3 (<i) 0-7 (^2) 



^AaA^CTa (A)rf / rf^l rft2 rfia /a (Xl,tl) //3 (Xi,t2) //3 (xi,t3)crQ (ti)o-/3 (i2)o-/3 (^3) , (59) 

R Jo 



r 



where tap-y is the antisymmetric tensor ll46ll . There is only 
one spatial index in ( l59l l because of the hypercube assump- 
tion: we have included only terms in which contraction of 
the /'s yields a non-zero value, as these will contribute to the 
effective short time evolution. At long times, connections be- 
tween the qubits are, of course, essential, and this is treated in 
the next section. 

The RG is naturally implemented in the case of ohmic baths 
(which leads to logarithmic singularities). However, suitable 
generalizations can be defined by dimensional regularization 
or by summing series in the expansion. Thus, in general, it is 
possible to write the following beta function for v'^^ : 



dXa 



■ 913-1 (^) XpX^i 



/3 



(60) 



where g and h are functions specific to a particular environ- 
ment, £ — A/A', and A' is the reduced (i.e. rescaled) cutoff 
frequency. By integrating the beta function from the bare cut- 
off, A, to A ~ ^ , we are summing the most relevant components 
of the noise inside a hypercube. If the renormalized value of 
the running coupling at frequency A^^, A*, is still a small 
number, then it is a good approximation to consider 



Va (Xi, A*) 



2h 



dt fa (Xi , t) 



(61) 



as 



the evolution operator of the qubit at position xi which was 



diagnosed with an error a by the QEC procedure. 

We illustrate the renormalization group procedure with two 
simple examples of ohmic baths: (i) marginally relevant and 
(ii) marginally irrelevant couplings. 

1. The k-channel Kondo problem 

The first example is a qubit exposed to a bosonic bath that 
is modeled by a SU{2)k Kac-Moody algebra - the bosonized 
Hamiltonian of a fc-channel Kondo problem. Here we closely 
follow the work of Affleck and Ludwig (see appendix B of 
Ref.liTi). We define chiral bosonic currents : Jl : obeying the 
operator product expansion (OPE) 



:J£(0::Ji(O: 



(62) 



where /'^''^ are the group structure constants and v is the ve- 
locity of excitations. In the interaction picture, the qubit cou- 
ples to the currents by the usual Kondo interaction, yielding 
an evolution operator (or, equivalently, a scattering matrix) of 
the form 



(63) 



Following our general discussion, we expand the evolution 
operator to lowest order in the coupling. 



U 



1-^/ dt:JL[t):-<j- 



a.h 



dt 



dt' ■.Jl{t)::Jl{t'): 



11 



(A \ 3 poo pt pt' 

-| ^ / dt dt' dt" : J£ {t) : : JI {t') : : J£ {t") : a VV^ (64) 

Due to the QEC evolution, only some of these terms are kept after the syndrome is extracted [see Eq. (I59ll1. For clarity, let us 
assume that we know from the syndrome that a phase flip has occurred. Hence, we must truncate the evolution operator to reflect 
this fact and apply the recovery operation (in this case multiply by a^), yielding 

■ \ pQO / \ \ 2 poo pt 

V. « I J^--J^L^^y--\2h) j Jtj Ji'\:-Jl{t)----Jl{t'):'':Jl{t)::Jlit')-] 

+ ^ (^) ' E j\ dt' dt" [: Jl it) : : J£ (f) : : J£ (t") : + : J^t) : : (f) : : J£ (t") : 

- ■.r^{t)::Jlit')::Jl{t"):]. (65) 



Now, we integrate over a small frequency shell [A — (5 A, A] 
and invoke the OPE. The result is a renormalization of the 
coupling A by an infinitesimal composed of quadratic and cu- 
bic terms. 



dX 



k . 



(66) 



The resulting running coupling X{£) can be used to improve 
the results of our bare perturbation theory. For that purpose, 
we integrate the beta function from the bare cutoff until A^^. 
For the case of a small number of channels, we obtain a renor- 
malized coupling of the form 



A* 



A 



1 - Aln|AA|' 



(67) 



Although the RG flow goes toward the strong coupling limit, 
we do not integrate the beta function all the way to zero fre- 
quency. Thus, if the renormalized coupling A* is still a small 
parameter, it replaces A leading to the first-order renormalized 
evolution 



iX*v 



/OO 
dt : Jl (t) : ct" 
-OO 



(68) 



2. Quantum frustrated system 

Correlations are not necessarily malignant to the com- 
puter's behavior. This is illustrate d by our second example: 
a quantum frustrated environment l48l l49l [s^ . Consider the 
case of three independent Abehan ohmic baths coupled as in 
Eq. but with the OPE 



■■Jl{t)::Jl{t'): 



2v2 [t - t'Y 



(69) 



Following precisely the same methodology of the previous ex- 
ample, we obtain the beta function 



dX 

m 



1 



A^ 



(70) 



which leads to the renormalized coupling 



A* 



A 



^1 + 2A2 ln|AA| 



(71) 



A quantum frustrated system has the remarkable property of 
asymptotic freedom. Hence, even very large bare couplings 
flow towards a perturbative regime. The physical reason be- 
hind this is the lack of a pointer basis 1151 11 . thus effectively 
decoupling the qubit from its surroundings ['49']. This phe- 
nomena can also be understood as self-inflicted 7r-pulse de- 
coupling working at the cutoff frequency A 10, E2l • 

If the three coupling constants have different bare values, 
then the flow stops at some finite frequency since two of the 
couplings will flow to zero before the third. In other words, 
there will be a pointer basis. In a quantum computer protected 
by QEC, however, we are effectively stopping the flow at a 
finite frequency. Hence, the effect described in the previous 
paragraph is relevant even for large anisotropic couplings. 



B. Probability of a faulty path 

Now that we have obtained a reasonable approximation to 
the evolution operator at each QEC step, we can turn to the 
problem of evaluating how much protection QEC yields at 
long times. The simplest quantity to calculate is the proba- 
bility of finding a particular history of syndromes, Eq. ( |54l ). 
Using Eq. ( |6T1 ) and the known commutation relations of the 
fa operators, we in general can write that 



and define 



< A,0 



.(iVA, (iV- 1)A) 



(A,0), (72) 



Xl.X 



dti dt2 fl^ (X, , i 1 ) fa,^ (Xj 



t2). 

(73) 

We now can evoke Wick's theorem once again to separate 
the intra- and inter-hypercube contributions to the probabil- 
ity: The quantum average (T^) ~ (i^ol '^w'^w |<<3o} canbe 
written as a sum of all possible pair contractions. It is conve- 
nient to separate the sum into two distinct parts. 
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First, the sum of all pair contractions in the same hyper- 
cube gives the stochastic error probability of a qubit, that we 
defined in Eq. (|49] |. namely, 

e„ = (xi, a;) |v3o) 

^ ^ dtidt2(/I(x,<l)/o(x,i2)), (74) 

where we used again that for |x — y| > ^ and \ti — <2| < 
A, we have (/l(x, <i)/o(y, <2)) ~ 0. Note that when 
we calculated A* we already summed intra-hypercube pair 
contractions; however, these were contractions on the same 
Keldysh branch [see Fig. [3ja)] and therefore are related to 
the wave function amplitude. Equation (|74| | corresponds to 
pair contractions between two distinct Keldysh branches [see 
Fig. Ob)], hence it gives the probability of that evolution. 
With this two-step procedure, we sum up the most relevant 
contributions to the probability within a hypercube. 

Second, we sum contractions between hypercubes. For 
each possible syndrome outcome we define the operators 



that is well discussed in the QEC literature. We rewrote the 
rest of the series taking into account the unusual non-unitary 
driven dynamics of QEC. The only remaining issue is to eval- 
uate the stability of the perturbation expansion in the renor- 
malized coupHng A*. 

In Sec. II Al we discussed how the scaling dimension of an 
operator is important when studying a perturbative expan- 
sion. The same argument holds when evaluating the protec- 
tion yielded by QEC in a coiTelated environment. If the scal- 
ing dimension of is 5a, then dim Fa = 25 a (see Appendix 
iGb . Hence, the original criterion for the validity of the pertur- 
bative expansion in X, D + z — 5a < 0, becomes 



D 



-25a <0 



(79) 



once the expansion in A* is adopted. Note the factor of 2 in 
this equation caused by QEC. 

Whenever Eq. ( l79l ) is satisfied, the long-range correlations 
will produce small corrections to the stochastic error proba- 
bility. Below, we illustrate this point with an example. 



Fo (x,0) = 1 ~ ^" ^^^^/'^^ : \fa (x,0)r : (75) 



and 



F„(x,0) = ^ (^^^ :|/„(x,0)|':, (76) 

where : : stands for normal ordering with respect to the envi- 
ronment ground state (see Appendix IFJi. We use these opera- 
tors to express the remaining pair contractions of each hyper- 
cube in the probabilities, namely. 



t'o'(x,A,0) 



i^o(x,0) 



and 



vl (x,A,0) «e„ [1 + (x, 0)] , 



(77) 



(78) 



Equations ( TTTT i and ( |78] l are the final ingredients needed to 
evaluate the probability of a particular history of syndromes, 
Eq. ( |54] |. The remarkable aspect of these equations is that they 
provide a very elegant reorganization of the perturbation se- 
ries. They were tailored to separate the local contribution, e^, 
from the long-distance, long-time components of the noise. 
Fa- The high-frequency part gives rise to the stochastic noise 



Probability of a "flawless" evolution. 

Consider the case of a non-Markovian noise model with 
only one type of eiTor (phase flips, for instance). For sim- 
plicity, assume that no spatial correlations exist {D = 0). 
Hence, we can consider each qubit separately and do not have 
to woiTy about the spatial structure of the quantum computer. 
We also assume a two-point correlation function of the form 



(/(x,ii) /(y,t2)) 



To 



h - to 



2(5/2 



(80) 



where tq is a constant with the dimension of time. How do 
these long-range correlations change the probability of a flaw- 
less evolution of a qubit after ^ 1 QEC steps? To answer 
this question, we evaluate 



N-1 



« {1-e f {^o\l[Fo{^i,jA)\^o}- (81) 

3=0 

Assuming e, A* ^ 1, we can rewrite the probability as 



(i^ol Tf exp ■ 



[A* A/ (2ft)]' 



JVA 
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1 + 



[{X*A/2h)f (To/A)4V^jV2(i-2a/^) 
(l-e)2 2{l-2d/z){1^4S/z) 



(82) 



r 



where we have kept only the leading term. There are two 
simple limits: 

(i) If z < 2(5, the corrections become increasingly irrele- 
vant as N grows. The stochastic probability in the limit of 



and the correction due to 



large N is given by 7-" (To) 
correlations are small. 

(ii) The tipping point is z = 26. By summing the subset of 
dominant terms 



NA 



dti 



A j„ 



dt 



2j+l 



|/(<2.-i)r::|/(t2.)r 



we obtain 



^(To) 



-Nf. 



1 



1 - 



In TV 



(83) 
(84) 



This signals a problem with the perturbative expansion when 
N « exp {fijr^) ■ For times larger than Aexp {fijr^) , 
correlations substantially change the probabiUty. 



C. Residual decoherence 

In addition to the probability of a given syndrome history, 
we also identified the residual decoherence, Eq. (ISST l. as a fun- 
damental quantity to QEC. The reason is that the noise models 
that we consider do not satisfy the Lafflame-Knill condition 
for perfect error correction 1 12], as is the case for most physi- 
cally relevant decoherence mechanisms. Hence, it may not be 
safe to ignore these high-order events in the couphng A. 



J 



It is straightforward to develop a calculation for the den- 
sity matrix along the same lines used for the syndrome history 
probability. After separating the intra- and inter-hypercube 
contributions, the perturbative expansion is reorganized using 
the renormalized coupling A* . The result is exactly the same 
as for the case of the probability: If D + z — 26 < 0, the per- 
turbation theory in A* is stable and the analysis of the residual 
decoherence done with the corresponding stochastic model is 
a good approximation of the true quantum result. We revisit 
the example used in Sec. lIIIBl to make this point clear. 



Decoherence of a "flawless" evolution. 

For this example, we assume an environment that can only 
introduce phase flip errors in the computer. As we discussed 
in Sec. U for this error model we can use the simple 3-qubit 
code. However, unlike the calculation of the probability of a 
flawless evolution, we now make some assumptions about the 
spatial structure of the computer: We consider for simplicity 
that each logical qubit is composed of three adjacent physical 
qubits. The encoding and decoding are described in Fig.[T] 

Following Ref. 39, we write the evolution operator for a 
particular logical qubit in a QEC cycle as 

wo(0,xo) = -uo (xi,0) Wo (x2,0)uo (x3,0) . (85) 



By expanding Eq. (l85b in powers of A, we obtain 



Wo (0,xo) 



2h 



E 



A r-ti 

dti / dt2f {Xj,ti) f {Xj,t2) 







3 



dti 



dt-. 







A 



dt3/(Xl,il)/(x2,ti)/(x3,ti)Z, 



(86) 



where Z is the logical phase flip for that particular logical qubit. Note that the third order term keeps the logical qubit inside the 
logical Hilbert space [39] and therefore is not corrected by the QEC code. 

We choose to evaluate the most off-diagonal term of the reduced density matrix. 



{^0\ (V'oinLAr-infcilWoO'A,Xfc 



it 



(<Po| (^ol rij^o^ rifeli Wo (xfc, jA) IV^o) l^'o) 



(87) 



where '\ — || ... t) and i = |i ... j) denote the state of the physical qubits, x^ is labeling M logical qubits, and N is the total 
number of QEC steps. 

After integrating all the modes inside a hypercube, we define a renormalized coupling A* and a local error probability e. 
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Finally, we evoke again Wick's theorem to write 

I- A- NMe^ ^e^{^y E..y /J"'" dt, J^' dt, {: f {^,h) : : f (y, t,) :) 



Pfj (T^o) = (V'oli) (Tl^o 



1-A + NMe^ + 64 (^)^ _ dti (: (x, : : (y, ^2) :) + 

where yl is a number proportional to e and A*. Hence, for e, A* <C 1, this simplifies to llssll 



(88) 



1 - 2NMe^ - 2e* 



A*A 



dx / dy 



NA 



dh / dt2(:/'(x,ti): :/2(y,i2):) + 



(89) 



If we now recall the two-point correlation function of Eq. (|7]), 
it becomes clear that the corrections due to correlations are 
relevant when D + z > 2S. 



D. Relation to the work of Aharonov, Kitaev, and Preskill 

The study of correlated noise has been a central problem 
for quite some time. Among the most recent advances is a 
paper by Aharonov, Kitaev, and Preskill (AKP) [261 . Using a 
method completely different from ours, AKP proved that: For 
a computer where qubits are interacting through an instanta- 
neous interaction of the form / Ax^^, it is possible to prove 
resilience for A < Ac and D — 26 < 0. The key distinction 
between the work of AKP and ours is the instantaneous na- 
ture of their interaction. Hence, while in our work each qubit 
is inside a distinct hypercube, for AKP they are all contained 
in a single hypercube. There is however a trade-off. Since 
their interaction is instantaneous and perfect error correction 
is assumed, there is no propagation of errors in time through 
the gauge field of the envkonment. Hence, effectively, AKP 
are considering a model with z = 0. As a result, our Eq. ( |79] l 
holds in the case they analyzed as well. 



V. THRESHOLD THEOREM AS A QUANTUM PHASE 
TRANSITION 

The main result of fault-tolerant quantum computation is 
the threshold theorem. The theorem states that if a stochastic 
error probability e is smaller than a critical value Ec, then the 
introduction of an additional layer of concatenation improves 
the protection of the information. Hence, for a fixed e, it is 
possible to sustain a quantum computation for any desire time 
at the cost of some reasonable additional hardware overhead. 

Even though quantum computation and QEC are out-of- 
equilibrium problems, it is intuitive to talk about different 
phases in the computer-environment parameter space. Along 
this line of thought, each phase corresponds to a distinct 
steady state. A natural choice for an order parameter is that 
given by the entanglement among the qubits and the environ- 
ment. We summarize our thinking in Fig.|5] where we present 
a schematic phase diagram for a quantum computer running 
QEC. 



For stochastic noise models, such an idea was explored by 
Aharonov ifTsll . Following that work, we can separate the be- 
havior of the computer into two distinct regimes: 

(i) For e < ec, the computer components can maintain 
large entanglement through fault-tolerant procedures, which 
in turns means that the computer and the environment are 
weakly entangled. Hence, due to this large internal entangle- 
ment, the quantum computer departs from the classical com- 
puter model. We can formalize these remarks by remember- 
ing that QEC tries to keep the system in the "steady state" 
described by the reduced density matrix. In order to keep the 
notation simple, lets take the ideal computer state as a pure 
state, 

= (90) 

with — J2i Q^i (t) I*} expressed in terms of the computa- 
tional basis {\i)}. As consequence, it has a reduced entropy 
5 « 0. In this case, if we look at the full Hilbert space (that is, 
before tracing out the environment), we find the tensor state 

1^-) « IV-) ® If^environment) • (91) 

(ii) For e > ec, the computer components are weakly en- 
tangled and, therefore, can be efficiently simulated by a Tur- 
ing machine. In other words, the computer density matrix no 
longer represents a pure state, but rather a statistical mixture. 
Thus, the computer components are strongly entangled with 
the environment. This corresponds to a steady state with a 
large reduced entropy (in the limit of e ^ 1, 5 « In 2, with 
N the number of qubits). 

In such a description, we see that e plays a role analo- 
gous to an effective temperature lf54ll . Hence, the threshold 
theorem defines a phase transition from a high-temperature 
phase, where qubits are independent from each other, to a 
low-temperature phase, where quantum coherence and entan- 
glement are possible [18]. This also sheds new light on the 
role of periodic measurements in QEC: They can be seen as 
a refrigeration that extracts entropy from the computer (very 
much like the Schulman- Vazirani initialization procedure ||55I1 
or the transfer of entanglement to fresh ancillas |3l). If the 
entropy production in the computer is below a certain level, 
then the computer can be kept in its "low-temperature" phase. 

Our analysis of correlated noise also fits perfectly into this 
description. The dimension criterion provided by Eq. ( |79l ) is 
the hallmark of a quantum phase transition It561 . For D + z < 
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FIG. 5: (Color online) Phase diagram of a quantum computer running QEC. The parameter 5 is the scaling dimension of the environment 
operator, D is the dimensionality of the computer, and z is the dynamical exponent of the environment [see discussion preceding Eq. (|79H. In 
the red phase, qubits and environment are strongly entangled causing strong decoherence. In the light blue phase, QEC keeps the qubits and 
environment disentangled, making computation possible. 



2(5, V can only produce small corrections to the stochastic er- 
ror model. The steady state of the system is therefore given by 
Eq. ( |9T| l. There is a clear separation of scales and the thresh- 
old theorem holds as it is. Conversely, for Z) + z > 2(5, there is 
no clear separation of scales. The computer and the environ- 
ment become increasingly entangled and the system is driven 
towards a different steady state. Such a state is probably dis- 
tinct from the "high temperature" one and it is likely that it is 
characterized by a smaller residual entropy. 

This does not mean that for + z > 2(5 it would not be 
possible to perform quantum computation. It only means that 
the threshold theorem as we stated it does not hold. It is con- 
ceivable that some different derivation of the theorem exists 
in this case. In this sense, 13 + z = 2(5 defines what is usually 
referred to as the upper critical dimension of the model (see 
Appendix IBTi. Below the upper critical dimension, there can 
be substantial corrections to the steady state given by Eq. W\\ , 
but it may still be possible to prove resilience. The question 
that remains open is whether a lower critical dimension exists, 
namely, a criterion for V that would tell us when it is impos- 
sible to perform long-time quantum computation. 



VI. SUMMARY AND CONCLUSIONS 

Most previous discussions of QEC have used the quantum 
master equation and quantum dynamical semi-groups |.57ll . 
This is a very natural approach: The computer is the object 
of interest; hence, one starts the discussion by integrating out 
the environmental degrees of freedom. However, the price 
paid in this approach is that some simplification is needed in 
order to derive the quantum master equation f^sl W\. The 
usual assumption is the Born-Markov approximation OZSll . In 
that case, it is natural to define an error probability for a given 
qubit, and a discussion in terms of error models naturally fol- 
lows fll l29ll . The situation is much less clear when the Born- 
Markov approximation cannot be justified I28i l58tl . In this 
case, temporal and spatial correlations can build up and com- 



pletely destroy the notion of the probability of an error. 

A key characteristic of the discussion here is that we do not 
try to use a quantum master equation. Rather, we follow the 
approach put forward by Schwinger and Keldysh ll42ll43i l44ll 
to study out of equilibrium systems. The main conceptual dif- 
ference is that we trace the environmental degrees of freedom 
only at the very last step of the calculation. Hence, we can 
make the most of the unitary evolution of a quantum mechan- 
ical system. 

Following this "Schwinger-Keldysh" approach, we dis- 
cussed the evolution of a quantum computer operated with 
fast and slow gates. On the one hand, for fast gates the mi- 
croscopic Hamiltonian is the one relevant for the evolution of 
the computer, Eq. ( l23T l. On the other hand, for slow gates we 
demonstrated that a suitable effective Hamiltonian, Eq. ( |40l ), 
can be used to provide an upper bound for the discussion of 
decoherence. With this effective Hamiltonian, the notation 
can be unified, and both cases treated simultaneously. We de- 
rived two formal expressions that quantify the evolution of 
the computer under QEC in a correlated environment: (i) the 
probability of a given syndrome history, Eq. ( l54b . and (ii) the 
reduced density matrix of the computer, Eq. ( |55] ). 

In order to fully use standard QEC theory, we introduced 
the important assumption of "hypercubes", that is a minimum 
spatial distance between qubits, Eq. ( fSTl l. in order to allow the 
definition of an error probability for a single qubit. With this 
"hypercube assumption", it is straightforward to use Wick's 
theorem to separate the environmental modes into intra- and 
inter-hypercube parts. The intra-hypercube component de- 
fines the error probability, while the inter-hypercube part is 
tracked by an operator acting on the coarse-grained scale of 
the hypercubes. As examples, we treated a generalization of 
the spin-boson model and a quantum frustrated model. 

All the pieces are put together when we explicitly calculate 
the probability of a syndrome history (Sec. lIVBl ) and associ- 
ated residual decoherence (Sec. lIVCl ). The main result is cast 
as a dimensional criterion, Eq. ( |79] ). Finally, we discuss the 
parallels between the threshold theorem and a quantum phase 
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transition. A qualitative description of the possible fates of a 
quantum computer as a function of noise strength and degree 
of correlation is given in Fig.|5] 

There are several clear directions in which our results could 
be extended or improved. First, it would obviously be desir- 
able to relax the hypercube assumption introduced in Sec. |lVl 
There is nothing intrinsic to our approach which makes this 
assumption necessary. Yet, progress without it seems much 
more difficult: The notion of a local error probability during a 
single QEC cycle becomes problematic, making the connec- 
tion with analysis based on error models, such as the usual 
derivation of the threshold theorem, unclear. 

Second, non-instantaneous gate operation is clearly a deli- 
cate issue. By using a bound (Sec. IllCb . we are able to treat 
this case in the same way as the fast-gate case. Thus we derive 
an upper bound for the local error probability together with the 
dimensional criterion. If a more accurate value for the error 
probability is desired, a specific error correction code as well 
as the gates under consideration must be included in the analy- 
sis. However, the scaling argument and resulting dimensional 
criterion do not, in general, change. 
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APPENDIX A: ABSOLUTE CONVERGENCE OF 
DYSON'S SERIES 

Dyson's series is absolutely convergent for any bound op- 
erator evolving for any finite time [62]. This is particularly 
simple to see using the sup operator norm ll28ll . 



\\A\\^snpJ{^\A^A\^), (Al) 

where = 1. If P = dt' \\V {t')\\ < oo, then the 

norm of the m*'^-order term in Dyson's series is bounded by 
P™ /ml. Thus, using the convergence of the exponential se- 
ries, we find that Dyson's series is absolutely convergent. 



Note that it is possible to change the dimensional criterion 
for the better (but not for the worse) by using the separation 
of scales introduced by QEC. Particular pulse sequences can 
reduce correlation at long times at the cost of increasing the 
local error probability. One example was given in our previous 

work mm- 

Finally, there may be a regime of parameters where, as indi- 
cated in Fig.|5] fault-tolerant quantum computation is possible 
even though the presently known derivations of the thresh- 
old theorem do not apply. By analogy with phase transi- 
tion phenomenology, there may be a lower critical dimen- 
sion such that a more sophisticated analysis than the one we 
present here shows that fault-tolerant computation is possible 
for d < {D + z)/2. It would be very interesting to show in any 
example that such is, or is not, the case. 

Quantum Error Correction is one of the most interesting 
frameworks which allows long quantum computations lf59tl . 
Even though QEC is widely accepted, it has bee n argued that 
it relies on a set of unphysical assumptions l25l l60l l6lll . 
namely: (i) "fast" measurements, (ii) "fast" gates, and (iii) de- 
scribing decoherence by error models. Although these are le- 
gitimate concerns, it is now clear that they are not fundamen- 
tal: First, in Ref. 35 DiVincenzo and Aliferis demonstrated 
that resilient circuits can be constructed with slow measure- 
ments. Second, in the current paper, we have demonstrated 
that the fast gate assumption is not critical for fault tolerance. 
Finally, we have laid the groundwork here for a theoretical 
framework that connects microscopic Hamiltonians with er- 
ror models in correlated environments. From our results for 
the threshold theorem in conjunction with those of AKP | |26l] . 
it is clear that a large class of correlated environments are al- 
ready properly treated within the QEC framework. 



APPENDIX B: PERTURB ATI VE EXPANSION IN 
(A* THEORY 



A classic example of a quantum phase transition is given 
by the 0"* theory at criticality [63i1 . The model is compactly 
described by the Euclidean action 



S = 



/3 



dT 



(Bl) 



The scaling dimension of the free field is usually defined as 
dim [0] = v/2. If we expand the partition function in pow- 
ers of A, it is simple to see that each order in the pertur- 
bative expansion will have the power \{L[3)^^^~'^ . Hence, 
D + l — z^<Ois the criterion for the irrelevance of the pertur- 
bation. The simplest way to see that is to do power counting 
by rescaling space and time. 



r — > 6r, T — > 5r, 
which immediately gives 



S = 6^-1-^ / d^r / dT 



b 2 , 



(B2) 



Xb 



D+1-2V / ,_D 



d^r / dT 



(B3) 



One finds the scaling A Afo^^^"", which is valid at each 
order of the perturbative expansion. The criterion for the ir- 
relevance of the perturbation is, D + \ — v <Q 

There is one more important definition that this example 
provides. Since the Gaussian action must be scale invariant, 
we automatically see that for this example v — D — 1. Hence, 
the criteria for the irrelevance of A^"* term as a perturbation 
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can be rewritten as 3 — D < 0. This defines the upper critical 
dimension for the model as d"PP°'' = 4 (three spatial and one 
temporal). When a system is above its upper critical dimen- 
sion, the physics is controlled by the Gaussian action. How- 
ever, when the system is below its upper critical dimension, 
there are substantial corrections to physical quantities when 
compared with the Gaussian solution. 



APPENDIX C: HILBERT SPACE OF QUBITS 

Due to the state vector normalization, the Hilbert space of 
a qubit is isomorphic to a three-dimensional sphere S^: For a 
general state = a |1) + /3 |0), we have the constraint 



{Reaf + (Imaf + [Re P) + (Im/3)^ = 1. 



(CI) 



However, an overall phase is physically irrelevant and the cor- 
rect mapping is to the complex projective plane of complex 
dimension 1, 



S^/U{1) CI 



(C2) 



For the same reason, the Hilbert space of n qubits is isomor- 



phic to C] 



])27l — 1 



For the discussion of entanglement, there 



is a particularly important subspace of this space. It is com- 
posed by the direct product of each qubit Hilbert space minus 
an over all phase. 



C 



(C3) 



modulus phase 

where j labels the j*'* qubit's Hilbert space. The dimension of 
the subspace grows as n — 1 while the dimension of the entire 
Hilbert space grows as 2n — 1. Entangled states are defined as 
the complementary set of this special subspace. 



APPENDIX D: DECOHERENCE IN THE SPIN-BOSON 
MODEL WITH OHMIC DISSIPATION 

An example of a qubit coupled to an environment is the 
spin-boson model with ohmic dissipation ll64ll65ll . which was 
intensively studied in the context of quantum computation 
l66ll even before quantum error correction was introduced. In 
this model, a qubit evolves according to the Hamiltonian 



H : 



dx {dx<i>Y 



A 9,0(0) a' 



(Dl) 



where </) is a chiral bosonic field, a are Pauli matrices that 
describe the qubit located at a; = 0, and A is the environment- 
qubit coupling constant. If a qubit is prepared in an initial 
state 



IV') = «|T) + /?|i), 



(D2) 



at large enough times, A ^ -C t <C (kBT) ^, its density 
matrix evolves as 



P{t) = 



_a*/3e-^'Mi+A*) |/3|2 



(D3) 



with A denoting the environment ultraviolet cutoff frequency. 
Since states with either a or f3 equal to zero do not experi- 
ence decoherence, they are called classical states. They de- 
fine a pointer basis. Conversely, any superposition state with 
a, /3 7^ suffers decoherence and over a long time becomes a 
statistical mixture of the classical states. 

As one includes more qubits, the entries in the reduced den- 
sity matrix will decay faster as one moves away from the di- 
agonal. In the case where qubits are coupled to independent 
baths, it is simple to see that the off-diagonal matrix elements 
decay as 



(t > A" 



-A^(p-9)ln(l+At) 



(D4) 



where p and oare the total magnetization of the states p and q, 
respectively [2]. The case of a common bath is also straight- 
forward, [2], and the result for qubits separated by a distance 
smaller than A^^ is 



-X^ip-qf ln(l+At) 



(D5) 



Some entangled states do not suffer decoherence (a sin- 
glet state, for example). However, these correspond to a 
very special and small decoherence-free subspace. In gen- 
eral, entangled states are made of quantum superpositions and 
therefore have components in the off-diagonal entries of the 
density matrix. Hence, studying decoherence (the decay of 
the off-diagonal elements of the density matrix) is essentially 
equivalent to studying how entanglement between qubits is 
destroyed by interaction with the environment. 



APPENDIX E: INTERACTION PICTURE 

Since [Hq, Hqc] = 0, we can define the interaction picture 

0{t) = ey^"°*R''{t)OR{t)e'-k"'>\ (El) 
= e-fe^°*i?t(t)f7(i) |*(0)), (E2) 

where U{t) is the exact evolution operator, defined as 

and I^I^) is the total state vector (computer plus environment). 
Now, let us consider the time evolution of |^), 



:V{t)\^it)) 



(E4) 



Thus, we obtain the usual definition for the evolution operator 
in the interaction picture 

Uit)=ei"°'R^{t)Uit)=Tte-i^o'^''^^''\ (E5) 
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APPENDIX F: LOW FREQUENCY CONTRIBUTION TO 
THE ERROR PROBABILITY 



The simplest way to understand F„ is to write / in its fre- 
quency representation 




^ K [ dt [/> (xi ,t) + f< (xi , 0)] , (Fl) 

where < stands for frequencies smaller than and > 
for the frequencies between and A. Thus, using that 

I 



(/a/a) =0. we obtain 

Jo 

+ (A:A)V<Uxi,0)/<(xi,0). (F2) 



APPENDIX G: SCALING DIMENSION OF F„ 

If the two-point correlation function of fa can be expressed 

as 

(/„ (xi,ii) /„ (X2,t2)) ~ ^ ((A^' Ja^) ' 

(Gl) 

the scaling dimension of /„ is defined as dim = S. Using 
Wick's theorem, 



(: |/a(xi,ii)|': :|/„(x2,t2)|' :) = (/l (xi, ti) /I (x2, t2)) (/« (xi, ii) (x2, is)) 

+ {fl (Xi, ii) fa (X2, t2)) (/a (Xi, ii) fl (X2, ti)) 

Therefore, dim Fa = 25. 
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